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: ABSTRACT 



We have monitored a sample of 27 nearby globular clusters in the 
hard X-ray band (20 - 120 keV) for ~1400 days using the BATSE in- 
> ; strument on board the Compton Gamma-Ray Observatory. Globular 

clusters may contain a large number of compact objects (e.g., pulsars 
or X-ray binaries containing neutron stars) which can produce hard X- 
ray emission. Our search provides a sensitive (~50 mCrab) monitor for 
0> ■ hard X-ray transient events on time scales of ;^ 1 day and a means for 
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observing persistent hard X-ray emission. We have discovered no tran- 
sient events from any of the clusters and no persistent emission. Our 
observations include a sensitive search of four nearby clusters containing 
c3 . dim X-ray sources: 47 Tucanae, NGC 5139, NGC 6397, and NGC 6752. 

The non-detection in these clusters implies a lower limit for the recur- 
^ ! rence time of transients of 2 to 6 years for events with luminosities 

^q36 gj,g g-i — 120 keV) and ~ 20 years if the sources in these 
clusters are taken collectively. This suggests that the dim X-ray sources 
in these clusters are not transients similar to Aql X-1. We also place up- 
per limits on the persistent emission in the range (2 — 10) x 10"^^ erg s~^ 
(2(7, 20 — 120 keV) for these four clusters. For 47 Tuc the upper limit is 
more sensitive than previous measurements by a factor of 3. We find a 
model dependent upper limit of 19 isolated millisecond pulsars (MSPs) 
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producing gamma-rays in 47 Tuc, compared to the 11 observed radio 
MSPs in this cluster. 

Subject headings: clusters:globular, 47 Tucanae, NGC 5139, NGC 6752, 
NGC 6397, X-ray:transients 



2 



1. Introduction 

Globular clusters (GCs) are remarkable 
stellar systems where a variety of compact 
objects may form and evolve. Numerous mil- 
lisecond pulsars (MSPs) inhabit GCs as re- 
vealed by radio observations (Lyne 1996), e.g. 
11 MSPs in 47 Tuc (Robinson et al. 1995). 
In the X-ray band, the number of bright 
sources (1 — 10 keV) per unit stellar mass 
is two orders of magnitude larger for GCs 
than for the remainder of the galaxy (Clark 
1975). Persistent X-ray sources in GCs di- 
vide into two luminosity groups: one with 
low luminosities (L^ ^ 10^^-^ erg s~^, 0.5 — 
4.5 keV), and the other with high luminosities 
(L^ 10^^ erg s-i) (Hertz and Grindlay 1983, 
Hertz and Wood 1985, Verbunt et al. 1995). 
The nature of the low luminosity, dim X-ray 
sources (DXSs) has been controversial (for re- 
views see Bailyn 1995, Verbunt et al. 1994, 
Grindlay 1993b, Hut et al. 1992). Sugges- 
tions for their origin include cataclysmic vari- 
ables, with mass accreting onto a white dwarf 
(Grindlay 1993a), low mass X-ray binaries 
(LMXBs) with accretion at quiescent levels 
onto a neutron star (NS) or black hole (BH) 
(Verbunt, Van Paradijs and Elson 1984), ple- 
rionic MSP binaries (Tavani 1991) or chance 
superpositions of foreground/background ob- 
jects (Margon and Bolte 1987). 

We have employed the BATSE instrument 
on board the Compton Gamma-Ray Obser- 
vatory to search for transient and persistent 
hard X-ray emission (;^ 20 keV). Hard X- 
rays may be produced by GC compact ob- 
jects in a number of ways. In a binary, hard 
X-ray emission can be powered by mass ac- 
cretion onto the BH/NS primary, or in pul- 
sars, from shock interactions in the relativis- 
tic pulsar winds in the binary. Hard X-ray 
emission may also be produced by magneto- 



spheric emission of isolated MSPs. Transient 
emission has been observed previously from 
GCs. NGC 6440, a GC containing DXSs, has 
exhibited a high luminosity transient episode 
(Forman, Jones and Tananbaum 1976). M15 
may also have shown a transient event (Pye 
and McHardy 1983). In the hard X-ray band, 
SIGMA has observed transient emission from 
the globular cluster Terzan 2. Recent ROSAT 
HRI observations indicate that the SIGMA 
source is most hkely associated with the X-ray 
burster XB 1724-30 (Mereghetti et al. 1995). 
SLX 1732-304 in the GC Terzan 1 has also 
been detected as a hard X-ray source (Chura- 
zov et al. 1994). Our search does not include 
Terzan 2, Terzan 1 or NGC 6440 due to prob- 
lems with nearby interfering sources. We have 
been able to observe M15 however. 

In Section 2, we describe hard X-ray emis- 
sion mechanisms in detail. Section 3 outlines 
the BATSE earth occultation flux measure- 
ment technique and the parameters of our 
search. In Sections 4 and 5 we present the re- 
sults of our search, providing upper limits on 
transient events and constraints on the recur- 
rence times of transients from DXSs. Section 
6 contains upper limits on persistent emission 
and constraints on the number of isolated and 
interacting MSPs. Section 7 is a summary 
and conclusion. 

2. Hard X-ray emission from compact 
objects 

2.1. Accreting black holes 

Soft X-ray transients (SXTs) have exhib- 
ited emission in excess of 10 keV (e.g. White, 
Kaluzienski and Swank 1984). These hard 
power law tails along with an anticorrelation 
of intensity and spectral hardness have been 
taken as indications of a BH primary (Tanaka 
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1989, Tanaka and Lewin 1995). 

Originally, the discovery of X-ray sources 
in GCs was interpreted as a possible mani- 
festation of massive accreting BHs (e.g., Bah- 
call and Ostriker 1975). We know today that 
all bright X-ray sources in GCs show X-ray 
bursts, and therefore most likely contain neu- 
tron stars. Though it is unlikely that clus- 
ters contain very massive black holes, they 
may harbor a population of stellar mass black 
holes. Kulkarni, Hut and McMillan (1993) 
estimate that there may be ~ 10 systems in 
galactic GCs in which a BH of mass ^ IOMq 
has captured a companion. In a few of these 
systems the companion may have had time 
to evolve off the main sequence, producing 
an active X-ray binary. Hard X-ray activity 
from such systems would be observable in our 
search. 

2.2. Accreting neutron stars 

LMXBs that contain NS primaries and ac- 
crete with luminosities below some critical 
value also exhibit transient episodes of hard 
X-ray emission (e.g. Barret and Vedrenne 
1994) and an anticorrelation of spectral hard- 
ness and intensity, e.g.: 4U0614-I-091 (Barret 
and Grindlay 1995a), 4U1608-522 (Mitsuda et 
al. 1989). 

Table |l| lists the observed hard X-ray prop- 
erties of several transient and persistent LMXBs. 
All of these sources have been observed as 
type-I X-ray bursters which indicate NS pri- 
maries. The hard X-ray outbursts generally 
exceed 10^^ erg s^^, with emission extending 
above 100 keV. The time scales for the hard 
X-ray outbursts are typically 10 — 100 days. 
Such events would be easily detectable in our 
search. Additional bursters are being de- 
tected as hard X-ray sources in a BATSE 
monitoring program (Barret et al. 1996). 



2.3. Isolated and binary Pulsars 

Persistent hard X-ray emission may also be 
produced by MSPs, and the fact that GCs 
contain large numbers of MSPs makes them 
ideal systems in which to search for such emis- 
sion (Chen 1991, Tavani 1993b). Recent the- 
oretical work has shown that hard X-ray and 
gamma ray emission may be produced by 
MSPs either in isolation (Chen, Middleditch 
and Ruderman 1993) or in interacting bina- 
ries (Tavani 1993a). Isolated MSPs could pro- 
duce hard X-ray emission by the same mech- 
anism operating in young Crab-like pulsars. 
Though the magnetic fields are smaller in re- 
cycled pulsars, the voltage drop at the light 
cylinder may be similar due to the smaller ro- 
tation periods (Ruderman and Cheng 1988). 
In binary systems, MSPs may emit hard X- 
rays by interaction with a companion (Tavani 
1991). This is in analogy to PSR B1957+20 
in which the pulsar is evaporating its compan- 
ion (e.g. Kluzniak et al. 1988). In such sys- 
tems the pulsar is partly or fully buried in a 
gaseous envelope originating from its compan- 
ion. The interaction of the MSP wind with 
the surrounding material produces energetic 
emission by synchrotron radiation or inverse 
Compton scattering of shock accelerated par- 
ticles. Emission, which can be explained by 
this model, has recently been observed in the 
1 — 200 keV range from the binary pulsar 
PSR 1259-63 near periastron (Grove et al. 
1995, Tavani et al. 1996). 

3. Cluster Search with BATSE 

The BATSE instrument was designed to 
provide continuous coverage of the entire sky 
in the hard X-ray /gamma-ray band from 20 keV 
— 2 MeV (Harmon et al. 1992), and has an 
optimal sensitivity in the range 20 — 120 keV 
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Source 


D 


Lx 


Observation 




(kpc) 


(erg s^"*^) 




Aql X-1 


2.5 


1 X lO''^^ 


BATSE (20 - 120 keV)'^ 


Cen X-4 


1.2 


7 X 10^6 


Signe 2MP (13 - 163 keV) ^ 


XB 1724-308 (Ter 2) 


14.0 


2 X 10^^ 


SIGMA (35 - 200 keV) ^ 


KS1731-260 


8.5 


2 X 10^^ 


SIGMA 


4U1608-522 


3.6 


6 X 10^^ 


BATSE*^ 


Al 742-294 


8.5 


6 X 10^6 


SIGMA' 



Table 1: Hard X-ray outburst luminosities of five X-ray bursters. is the peak outburst lumi- 
nosity for the 20 — 120 keV band. The energy bands of each instrument are shown. 



"(Harmon et al. 1996) 
''(Bouchacourt et al. 1984) 
(Barret et al. 1991) 

''(Barret et al. 1992) 
^ (Zhang et al. 1996) 
^ (Churazov et al. 1995) 



for a source with a typical photon index of 2. 
Source fluxes are measured by BATSE using 
the earth as an occulting object. As a source 
rises or sets, steps appears in the detector 
count rates due to attenuation by the Earth's 
atmosphere. Source fluxes are determined by 
fitting the amplitude of these steps. Thus, de- 
pending on the geometry of the Earth's rising 
and setting limbs, one can make two source 
measurements per GRO orbit (approximately 
90 min.), and from these construct a light 
curve. 

Fluxes in the resulting light curve, for a 
location in the sky with no sources, form a 
distribution about a mean flux of zero. Both 
positive and negative fluxes are produced in 
the occultation edge fits, all with formal er- 
ror estimates. The distribution of flux esti- 
mates has a width which is larger than ex- 
pected for a Poisson distribution of back- 
ground rates. This is due to systematic ef- 
fects, such as contributions in the occultation 



flts from nearby sources and spurious back- 
ground fluctuations. 

The main limitation of the occultation tech- 
nique is source confusion due to the large an- 
gular extent of the Earth's occulting limbs. 
There are several diagnostics that can be em- 
ployed to identify confused sources. One such 
check is to compare the measured relative 
rates in separate detectors to the rates ex- 
pected given the BATSE detector responses 
(Pendleton et al. 1996). Another technique 
is to compare the rates from the rising limb to 
those from the setting limb. For constant flux 
from the expected source, the rising and set- 
ting edges should give the same rates. The 
precession of the earth's limb can also be 
used to localize the emission, a technique that 
has been successfully employed to create im- 
ages of sky regions for brighter, long duration 
events (Zhang et al. 1993). 

We divide our search of 27 GCs into two 
parts, which we will refer to as Search A 
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and Search B (see Figure |T^). Search A is a 
lower sensitivity search of 23 GCs which are at 
larger distances and/or in regions of the sky 
prone to source confusion. Search A includes 
GCs at distances ^ 6 kpc, GCs near the galac- 
tic center, GCs containing no known DXSs, 
and GCs known to contain bright LMXBs. 
Search B is a higher sensitivity search of four 
nearby GCs which are not subject to seri- 
ous source confusion problems. We limit this 
search to GCs at ^ 6 kpc which also contain 
DXSs near the cluster core. We also exclude 
clusters within 15° of the galactic center, to 
avoid source confusion in this dense area. The 



GCs of Search B are shown in Figure lb 



cir- 



cled). In what follows, all cluster parameters 
have been taken from Peterson (1993). 

4. Results from Search A 

We have generated light curves for all the 
clusters extending from April 1991 to March 1995 
(~1400 days). For Search A, we have an- 
alyzed the rate history of each cluster for 
events with structure on a time scale of one 
or more days. Generally in each light curve 
there are several "events" apparent; most of 
these can be attributed to interference from 
nearby bright sources. There are, however, a 
few spikes that cannot readily be explained as 
interfering sources. None of these events are 
particularly outstanding in terms of ampli- 
tude or structure. For Search B we have fur- 
ther analyzed such events to determine if they 
originate from the clusters, as described be- 
low. For the 23 clusters in Search A, however, 
such an analysis is impractical. For Search A, 
therefore, we take the largest of the outstand- 
ing features to set conservative upper limits 
to the minimum observable event amplitude 
(Table 0). We initially find count rate upper 
limits {Rui in Table To convert these lim- 



its to photon fluxes, we have multiplied by 
a constant conversion factor determined from 
other light curves for which we have properly 
deconvolved the instrument response with an 
assumed photon index of 2.0 to generate pho- 
ton flux histories. This conversion is accurate 
to within approximately 30%. From the pho- 
ton rate, we calculate a luminosity upper limit 
by assuming power law spectra with index 2. 
The effective time coverage for each cluster is 
shown in Table |[ Coverage with BATSE is 
less than 100% and varies as a function of sky 
position due to the geometry of the Earth's 
occulting limbs and due to interfering sources. 
When target occultations occur within 10 sec- 
onds of the occultation of a bright interfering 
source, no edge fits are produced. This leads 
to data gaps, as in the case of NGC 6760, 
which has highly non-uniform coverage due 
to the nearby source Aql X-1. Note that we 
obtain a limit of ~ lO^^ erg s'^ for NGC 5272 
(M3) which is known to contain a time vari- 
able supersoft X-ray source. Note also that 
the cluster NGC 6440 has been excluded from 
our search due to nearby interfering sources. 

5. Results from Search B 

For Search B, we have produced rate histo- 
ries over the time period mentioned above and 
properly deconvolved the detector response to 
obtain photon light curves. We have searched 
these light curves for outburst events. Can- 
didate events are identified as features of du- 
ration one or more days which rise at least 
1.5(7 above the background flux distribution. 
With this criterion we have identified 40 pos- 
sible candidate events. We analyzed each of 
these isolated events further to determine if 
they originate from the clusters. We check 
for nearby interfering sources and compare 
the relative rates in the detectors. This anal- 
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Cluster 


JJ 






T 


coverage 


Persistent 




(kpc) 


(cnt s ) 


(7 cm s ) 


(erg s^^) 


(%) 


X-ray 








\ 1 n— 21 
[xiU J 


Fx^ 1 n361 
[X iO J 




Source(s) 


AAA 

iV14 


2.U 


3.5 


2.8 


A A 

0.9 


OA A 

80.0 




IN (jUd544 


2.5 


4.0 


3.2 


1.7 


77.1 




iNLrCDDOD 


3.U 


2.6 


2.1 


1.6 


OO A 

82.4 


dim 




3.9 


/I A 

4.0 


3.2 


A A 

4.0 


OA r 

80.5 




NGC6539 


4.0 


A A 

4.4 


3.5 


A ''7 

4.7 


'^A A 

70.9 




iNlLrUDdDD 


4.U 


6.5 


5.2 


C A 

6.9 


'1 A A 

74.9 




NGrC6760 


4.2 


6.3 


r A 

5.0 


7.2 


o o ^ 

38.7 






A O 

4.3 


3.3 


2.6 


y< A 

4.0 


^ A O 

74.3 




i\IGGo8U9 


A O 

4.8 


4.5 


3.6 


6.8 


72.3 






5.U 


T A 

7.0 


5.6 


11.6 


77.0 




Aipp/l Q79 


o.z 




O.Z 


7 9 


OtJ.o 




NGC6218 


5.6 


3.8 


3.0 


7.8 


74.5 




NGC4833 


5.8 


2.9 


2.3 


6.4 


69.0 




NGC6626 


5.9 


4.0 


3.2 


9.2 


79.4 


dim 


NGC6541 


6.6 


3.3 


2.6 


9.3 


73.3 


dim 


NGC6712 


6.8 


4.4 


3.5 


13.4 


72.9 


bright 


NGC7099 


7.4 


3.0 


2.4 


10.9 


88.4 


dim 


NGC6341 


7.5 


3.5 


2.8 


13.0 


85.7 


dim 


NGC6624 


8.1 


3.5 


2.8 


15.2 


84.2 


bright 


NGC5272 


10.1 


2.4 


1.9 


16.0 


88.0 


dim 


M15 


10.5 


3.5 


2.8 


25.5 


87.6 


bright 


NGC1851 


12.2 


4.5 


3.6 


44.3 


88.1 


bright 


NGC1904 


13.0 


3.1 


2.5 


34.9 


88.1 


dim 



Table 2: Event upper limits and coverage time for the 23 GCs in Search A. R^i is the raw count 
rate upper limit for hard X-ray transient events, and F^i the approximate corresponding limit in 
flux units. The luminosity upper limit (20 — 120 keV), L^i, assumes the stated distances and power 
law spectra with index 2. The total time span is 1630 days for NGC 6544 and NGC 6656 and 1490 
days for the other clusters. The clusters containing bright and dim X-ray sources are marked. 
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ysis eliminates all but 13 of the events as 
candidates. These 13 candidates consist of 
peak-like structures with durations of 1 to 8 
days. The significance of these events, as mea- 
sured from the photon flux error bars near 
the peaks, ranges from 2a to Qa. The max- 
imum luminosities in these candidate events 
would range from 0.7 — 6 x 10^^ erg s~^, and 
in one case 9 x 10^^ erg s~^ (20 - 120 keV). 
These luminosities are similar to the expected 
transient peak luminosities (c.f. Table |l|) and 
deserve further investigation. 

We further analyze candidate events by 
studying the correlation between the rising 
and setting rates within a given GRO pointing 



period. Figures ^ and |b| show example rate 
histories for candidate events in NGC 6752 
and NGC 5139. For each source, the rise and 
set rates clearly have different temporal be- 
haviors. In each case nearly all the flux is 
from the fits in the rising limbs. Formally 
we can calculate the correlation coefficient for 
the rising vs setting rates. As a test of this 
technique we have calculated the correlation 
coefficient for a two day data sample from the 
Crab. The rise and set data for the Crab are 
well correlated with a correlation coefficient of 
0.77 and a corresponding probability of 2.6% 
that the observed correlation could arise from 
a random sample. For the two events shown 
in Figures |2^ and the correlation coeffi- 



cients are -0.07 and -0.10 respectively, with 
a 54% and 49% probability of being random. 
We have calculated the correlation coefficient 
for the other 11 candidate events during the 
relevant time intervals. For each event the 
rises and sets either show very weak correla- 
tion or, in several cases, strong anticorrela- 
tion. The most significant correlation has a 
30.3% probability of being random. The lack 
of significant positive correlations imply that 



the "events" observed in the light curves do 
not originate from the target sources but are 
due to interfering sources located elsewhere. 

As a check to this interpretation we have 
analyzed the event in NGC 6752 shown in 
Figure ^ in more depth. We have chosen 
to look further at this particular event be- 
cause it is somewhat extended in time, lasting 
for approximately 5 days, and reaches mod- 
erate fiux levels. We have produced a map in 
an approximately 7° x 7° region centered on 
NGC 6752 during the time of the candidate 
event. The map is created by forming a grid 
of points separated by 0.5°, and producing oc- 
cult ation histories for each grid point. The re- 
sulting rates for each point can be combined 
together to form a map. The map for the 
NGC 6752 event shows an emission enhance- 
ment well localized at ~ 3° 'northeast' of the 
cluster during the first GRO pointing period 
(TJD 9078 - 9083). It is also clear from the 
limb geometry that such an event would be 
visible only on the setting rates, and indeed 
this is what is observed in the rate history for 
NGC 6752. During the second pointing pe- 
riod (TJD 9083 - 9090), the emission peak 
moves to ~ 1° 'southwest' of the cluster, and 
is well localized only in the detector with the 
less optimal pointing. All of this leads us to 
conclude that the enhancement is not associ- 
ated with the cluster. 

5.1. Recurrence time limits from DXSs 

In the absence of any detected transient 
events, we are able to set a lower limit to the 
recurrence times of transient sources in the 
clusters of Search B. For a given time interval 
in a GC light curve, there will be a minimum 
luminosity, Lmm, at which an event would be 
visible. Lmin is determined by the width of 
the fiux distribution. We step through the 
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cluster light curves and determine the total 
number of days that events of different Lmin^ 
would be observable. The result is a distri- 
bution of the number of days observable vs 
Lmin (Figure |3^). The Lmin distribution de- 
pends on the assumed duration of outburst 
events and the photon index used in the de- 
convolution. This dependence, however, is 
rather weak. The distribution depends more 
strongly on the time intervals into which the 
data are integrated. Figure ^ shows the 
Lmin distribution for different data integra- 
tion times. The minimum observable lumi- 
nosities are shifted down at longer integration 
times due to an increased sensitivity, i.e. a 
smaller width in the distribution of fluxes. 

We determine limits to the recurrence times 
of outbursts from the Lmin distribution. For 
each Lmin there is an associated observable 
time, Tobs, for which an event exceeding Lmin 
would be detectable. Tobs is essentially the in- 
tegral of the Lmin distribution. Tobs is related 
to the mean event recurrence time, Tevent, by 
P = e~^°''»/'^'='"="*, where P is the probability of 
observing no events in time Tobs assuming a 
Poisson distribution of events. We calculate 
a lower limit for Tevent using a 5% probabil- 
ity of detection. With multiple X-ray sources 
in each cluster, an event could originate from 
any of the sources. The minimum event recur- 
rence time, Tevent, cau be stated as a minimum 
recurrence time for sources, if we use the num- 
ber of DXSs in each cluster, N^xs- The min- 
imum source recurrence time, r^, is given by 
Tr = Njjxs- Tevent (Table Wc also calculate 
a total recurrence time for all the sources as a 
class using the total Tobs for all clusters (which 
gives the total Tevent) and the total N^xs of 
15. Figure Ed is plot of vs Lmin for the 
four clusters. Recurrence times for sources in 
the clusters are constrained to lie above the 



lines. It is clear that with increasing event 
luminosity, Tobs increases and therefore in- 
creases, giving a tighter recurrence time con- 
straint at higher luminosities. Table ^ dis- 
plays recurrence time lower limits for event 
luminosities of 5 x 10^^ and 1 x 10^^ erg s~^. 
The estimates of N^xs are taken from the 
literature: 47 Tuc (Hasinger, Johnston and 
Verbunt 1994), NGC 5139 (wCen) (Cool et 
al. 1995a), NGC 6397 (Cool et al. 1993), 
and NGC 6752 (Grindlay 1993b). For 47 Tuc 
we take only the DXSs within the core. For 
the latter three clusters we have assumed, for 
consistency, that sources within 5 core radii 
of each cluster center are associated with the 
GC and others are not. Note that changing 
Ndxs simply scales the numbers in Table ^ 
and the curves in Figure |3c . 



These recurrence time limits can be com- 
pared to the several recurrent transients that 
have been discovered with all-sky X-ray sur- 
veys. One such transient is Aql X-1, observed 
over a 7 year period with the Vela 5B satel- 
lite. It has an outburst recurrence time mea- 
sured at 1.2 years (Priedhorsky and Terrell 
1984). Another example is Cen X-4 from 
which hard X-ray outbursts were observed in 
1969 and 1979 (Kaluzienski, Holt and Swank 
1980, Bouchacourt et al. 1984), indicating 
that its recurrence time is approximately 10 
years or less, putting it at the outer edge of 
our limit. Both Aql X-1 and Cen X-4 have ex- 
hibited hard X-ray emission during outburst 
in excess of 1 x 10^^ erg s"^ (20 - 120 keV) 
(c.f. Table |ip. A number of other sources have 
also displayed recurrent outbursts: 1608-522, 
1630-472, 1730-335 (Van Paradijs and Ver- 
bunt 1984, Chen, Shrader and Livio 1996). 
These systems have recurrence times of ap- 
proximately 0.5 to a few years. 

Our recurrence limits indicate that tran- 
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Cluster 


D 

(kpc) 


Ndxs 


(years) 

5 X 10^^ erg s"^ 


Tr 

(years) 


(years) 

1 X 10^^ erg s"^ 


(years) 


47 Tuc 


4.6 


5 


2.0 


3.4 


3.5 


5.8 


NGC 5139 


5.2 


2 


0.6 


0.4 


2.9 


1.9 


NGC 6397 


2.2 


5 


3.6 


6.0 


3.6 


6.0 


NGC 6752 


4.1 


3 


3.2 


3.2 


3.6 


3.6 


Total 






9.4 


13.0 


13.6 


17.4 



Table 3: Recurrence time limits for the GCs in Search I. Tobs are the integrated times per cluster 
during which events of the specified luminosity would be detectable. Tr are the corresponding max- 
imum recurrence time of outbursts from the dim X-ray sources (DXSs). N^xs are the estimated 
number of known DXSs in each cluster. 



sient sources with outburst luminosities greater 
than approximately 1 x 10^^ erg s~^ (20 — 
120 keV) and recurrence times less than about 
2 to 6 years cannot constitute the population 
of dim GC sources in the four sampled clus- 
ters. The recurrence time lower limit for all 
sources is ~ 20 years, if one is willing to con- 
sider them as a class. This suggests that the 
dim sources in these clusters are not quiescent 
LMXBs of a type similar to Aql X-1 or Cen X- 
4. Some caution should be taken with this in- 
terpretation. The outburst recurrence times 
are determined well only for a few systems 
which may not be typical. Also LMXBs in 
GCs may have systematically different prop- 
erties than the field binaries in which the hard 
X-ray emission and recurrent outbursts have 
been observed. Indeed some bright LMXBs 
in GCs have unusual properties, for exam- 
ple very short orbital periods (Bailyn 1996). 
One might argue that a substantial number 
of X-ray outbursts from Aql X-l-like LMXBs 
might be detectable only below 20 keV. If 
this were the case, such events would be not 
detectable by BATSE. However, the recently 
monitored behavior of Aql X-1 shows clearly 
that BATSE can usually detect hard X-ray 
emission at times when the optical counter- 



part of Aql X-1 is excited. These events are 
most likely associated with major X-ray out- 
bursts (Harmon et al. 1996). 

6. Limits from the Persistent Emission 
Search 

We have searched for longer time scale per- 
sistent emission using the flux histories for the 
four GCs in Search B. Previous upper limits 
have been obtained only for 47 Tuc. Obser- 
vations with SIGMA spanning seven days set 
a limit at ~ 2 x 10^^ erg s'^ (40 - 100 keV) 
(Barret et al. 1993) and an observation by the 
balloon-borne EXITE instrument found an 
upper limit of 1.4 x 10^*^ erg s"^ (30-150 keV) 
(Grindlay et al. 1993c). Upper limits have 
also been obtained for 47 Tuc with COMP- 
TEL, OSSE and EGRET in their respective 
bands (O'Flaherty et al. 1995). 

Our BATSE search offers improved sensi- 
tivity result of the very long integration 
time. We have detected no emission for the 
four clusters. The upper limits are summa- 
rized in Table ^. We have obtained these up- 
per limits by formal error calculations inte- 
grating over large time windows. As discussed 
in Section 3, such a procedure yields error val- 
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Cluster 


Ful 

(7 cm~^s^"'^) 
[xlO-^] 


Lui 

(erg s~^) 
[xl034] 


Nisol 


Finter 

(7 cm~^s^^) 
[xlO^^] 


■Winter 


47 Tuc 


3.3 


5.8 


19 


0.71 


46 


NGC 5139 


4.1 


9.2 


31 


0.62 


66 


NGC 6397 


4.2 


1.7 


6 


3.10 


14 


NGC 6752 


3.6 


5.0 


17 


0.85 


42 



Table 4: Persistent emission upper limits for GCs in Search I. The table lists Fui (values are 
~1 mCrab), the 2cr flux upper limits (20 — 120 keV), Luu the corresponding luminosity upper 
limits for power law spectra with index 2, Nisou the maximum number of isolated MSPs producing 
hard X-ray emission based on the model of Chen (1991), Finter-, the estimated flux produced by 
a single interacting MSP, and Ninter, the maximum number of interacting binary MSPs for an 
assumed pulsar spin-down power of 10^^ erg s~^ (Tavani 1993b). 



ues systematically larger than the actual flux 
distribution. We have studied this effect for 
the four GCs of Search B and determined a 
correction factor of 1.6 based on fitted widths 
of flux distributions. We have used the fac- 
tor of 1.6 correction to the formal error esti- 
mates for the upper limit values of Table ^ 
The clusters of Search B were chosen to limit 
the systematic effects of interfering sources, 
a selection which permits very sensitive flux 
limits. We note that the sensitivities may de- 
pend on the assumed spectral shape used in 
deconvolution. We have employed power law 
spectra with photon index 2.0. 

Using our upper limits, we constrain the 
number of isolated high energy emitting MSPs 
in the clusters. The luminosity for a typi- 
cal non-interacting MSP can be calculated. 
Chen (1991) does this using the outer-gap 
emission model of Cheng, Ho and Ruderman 
(1986) and the observed period distribution 
for MSPs. The resulting luminosity from a 
cluster is a constant times the number of iso- 
lated MSPs, assuming isotropic emission. For 
a cluster with one MSP, this would yield a lu- 
minosity of approximately 3 x 10^^ erg s~^. 



Combining this luminosity estimate with our 
luminosity upper limit, we find the maximum 
number of isolated MSPs, Nisoi, in each of 
the clusters (Table H). The estimate of Nisoi 
depends on the period distribution of MSPs. 
In the luminosity calculation, Chen (1991), 
uses the period distribution of all the 28 then- 
known MSPs, fit falling power law. This 
distribution is poorly constrained due to se- 
lection effects and small number statistics and 
may vary from cluster to cluster. This adds 
additional uncertainty to our values of Ni^oi- 
The flux from an interacting binary MSP 
has been estimated by Tavani (1993b). Us- 
ing this model we calculate the estimated 
flux from each cluster for a single interacting 
MSP system, F^ter (Table ^. Combined with 
our flux upper hmits, the flux estimate yields 
an upper limit on the number of interacting 
MSPs, Ninter- The constraints on N inter are 
weak. This calculation uses an average MSP 
spin-down power of 10^^ erg s~^ (Taylor et 
al. 1995). MSP spin-down powers are in the 
range lO^'^-lO'^^ erg s~^ but are still not well 
constrained. 
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7. Conclusions 

We have monitored 27 galactic globular 
clusters (GCs) with BATSE during a pe- 
riod of approximately four years each. We 
have detected no distinct hard X-ray outburst 
episodes and no persistent emission. The lack 
of detected events has several interesting im- 
plications. 

For the GCs with dim X-ray sources (DXSs) , 
we find a lower limit for the outburst recur- 
rence time from DXSs of ~ 2 — 6 years. 
This limit excludes the existence of 'Aql X- 
1-like' objects (i.e., persistent X-ray sources 
subject to major X-ray outbursts with a time 
scale of ~ 1 year), since is comparable to 
or greater than this outburst recurrence time 
scale. This suggests that the DXSs in these 
clusters are not LMXBs similar to Aql X-1. 

The lack of strong outburst events in our 
search means we also have no evidence of 
accreting BHs in GCs. If active BH bina- 
ries exist in GCs, they would be clearly de- 
tectable with hard X-ray luminosities of order 
of 10=^6 - 10=^^ erg s-^ 

We have calculated upper limits on persis- 
tent hard X-ray emission from GCs. The lim- 
iting luminosity is ~ (2 — 10) x 10^^ erg s~^ for 
the closest clusters. The limit for 47 Tuc is 
somewhat more sensitive than previous mea- 
surements. A model-dependent but reason- 
able estimate of the expected hard X-ray 
magnctosphcric emission from isolated MSPs 
(Chen 1991) implies that the number of iso- 
lated MSPs emitting hard X-rays in 47 Tuc is 
less than 19. Taking the observed 11 MSPs 
in 47 Tuc, and a beaming factor of 2 for short 
period pulsars, this upper limit is comparable 
to actual number of MSPs in 47 Tuc. Our 
results are therefore a constrain the magneto- 
spheric model of MSP hard X-ray emission. 



Our upper limits on persistent hard X- 

ray emission also provide a mild constraint 
on the number of interacting pulsars in bina- 
ries. The observed efficiency (a few percent) 
of conversion of pulsar spindown luminosity 
into hard X-ray emission in the case of the 
periastron passage of the Be star/pulsar sys- 
tem PSR B1259-63 (Tavani et al. 1996, Grove 
et al. 1995) together with the observed MSP 
spindown average luminosity of 10^^ erg 
implies a limit to the GC population interact- 
ing pulsar binaries. We find limits of 14 — 66 
interacting MSPs in the clusters studied. 

We would like to acknowledge the BATSE 
instrument team for their support. We thank 
Jonathan Grindlay, Didier Barret and An- 
drew Chen for helpful comments. We thank 
Mark Stollberg for assistance with data analy- 
sis. This work is Columbia Astrophysics Lab- 
oratory Contribution Number 596. This work 
is supported in by NASA Grants NAG 5-2235 
and NGT 8-52806. 
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Fig. la. — Projection of the positions of globular clusters onto the galactic plane. All the circled 
clusters are included in our Search A or Search B. 
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Fig. lb. — Globular clusters within 6 kpc projected onto the plane. The four clusters for which we 
have performed the sensitive Search B are circled. All the other clusters are included in Search A. 



16 





RISES 














\ 










t T - 

9086 - 9096i 

<R> = 1.12 ± 0.34 


- SETS 

\ 




\_ 


\- 


\- 






\\ 










9086 - 9096: 



9080 9090 9100 

TJD 



Fig. 2a. — Rate histories during candidate event in NGC 6752 near TJD 9090. The rate histories 
have been separated into rise and set rates, and corrected approximately for detector response. 
The average rise and set rates during the approximate time of the event are shown. 
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Fig. 2b. — Rate histories during candidate event in NGC 5139 near TJD 8990. 
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3a. — Coverage for the four clusters of Search A. Data integration time is 10.0 days. 
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Fig. 3b. — Coverage for the cluster 47 Tuc using various data integration times. 
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Fig. 3c. — Recurrence time lower limits for dim X-ray sources. Recurrence times are constrained 
to lie above the curves. A data integration time of 10.0 days has been used. 
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